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Abstract 

We study the exclusive decay Bg — > r+r^ in the general two Higgs doublet model. We 
analyse the dependencies of the branching ratio on the model parameters, including the 
leading order QCD corrections. We found that there is an enhancement in the branching 

ratio, especially for rth = > 1 case. Further, the neutral Higgs effects are detectable 
for large values of the parameter C^^^. 
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1 Introduction 



The study of rare B-decays is one of the most important research areas in particle physics and 
there is an experimental effort for studying them at various centers such as SLAC (BaBar), 
KEK (BELLE), B-Factories, DESY (HERA-B). In the Standard model (SM) they are induced 
by flavor changing neutral currents (FCNC) at loop level and therefore they are sensitive to the 
fundamental parameters, like Cabbibo-Kobayashi-Maskawa (CKM) matrix elements, leptonic 
decay constants, etc. These decays also provide a sensitive test to the new physics beyond the 
SM, such as two Higgs Doublet model (2HDM), Minimal Supersymmetric extension of the SM 
(MSSM) 0, etc. 

Among the rare B decays, l^l^ process, induced by the inclusive b — » sl^l^ decay, is 

attractive since the only non-perturbative quantity in the theoretical calculation is the decay 
constant of Bs which is reliably known. From the experimental point of view, the measurement 
of the hadronic decay is easier compared to its inclusive channel. 

The measurement of upper limit of Bg f^^fJ'" 

BR{Bs < 2.6 10^*^ , (1) 

stimulated the study of Bg l^l^ decays. In the literature, this process (0-0) and its 
inclusive one Bg —>■ Xgl^l^ (||TU|-||T^) have been investigated extensively in the SM, 2HDM and 
supersymmetric model (SUSY). When I = e, fi, the neutral Higgs boson (NHB) effects are safely 
neglected in the 2HDM because they enter in the expressions with the factor me{fj,)/mw, mH±. 
However, for / = r, this factor is not negligible and NHB effects can give important contribution. 



In [0, B XgT^T^ process was studied in the 2HDM and it was shown that NHB effects 
are sizable for large values of tan(3. Therefore the main observation of these calculations is the 
enhancement of the branching ratio {BR) of these decays for large tan(3 values in the 2HDM 
and minimal supersymmetric model (MSSM), especially for / = r lepton case. In a recent 
work [0, the inclusive h st^t~ process has been studied in the general 2HDM with real 



Yukawa couplings and it was found that the BR has been enhanced for large values of the 
parameters ^^^^ and ^^f,^. In this work, we study the Bg —>■ t^t^ decay in the general 2HDM, 



so-called model HI . Our calculations are based on the results of the work |]I4[ for the inclusive 
h ST^T~ decay. Here we include NHB effects and make the full calculation using the on-shell 
renormalization prescription. The investigation of the dependencies of the BR on the model 
parameters, namely ^^bf, and ^^^-t-, shows that a large enhancement in the BR is possible. 
The paper is organized as follows: In Section 2, we present the leading order (LO) QCD 
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corrected effective Hamiltonian and the corresponding matrix element for the exclusive Bg 
T^T~ decay in the framework of the model III. Section 3 is devoted to the analysis of the 
dependencies of the BR on the the Yukawa couplings ^^ (,6' ^n,tt discussion of our 

results. In Appendices, we present the operators appearing in the effective Hamiltonian and 
their Wilson coefficients. 

2 The Bs T^T~ decay in the framework of the model 
III 

The general Yukawa interaction in the model III is 

where L and R denote chiral projections L{R) — 1/2(1 ^ 75) and 0^ for i — 1,2, are two scalar 
doublets. Here rjfj^, are the Yukawa matrices and, in general, they have complex entries. 
The choice of scalar Higgs doublets 



\ ^x' 

(3) 



with the vacuum expectation values. 



<^^>=^[l) ■.<<t>2>=0, (4) 

and the gauge and CP invariant Higgs potential which spontaneously breaks SU{2) x U{1) 
down to U{1): 

V{<f>l,<f>2) = Ci(0f(/.i-^;72)2 + 02(02+02)' 

+ +C3[(0|0i - vy2) + 0+02]' + C4\{(pt(pi){(pth) - {(Pth){(pt(pi)] 

+ C5[i?e(0+02)]' + C6[/m(0+02)]' + C7, (5) 

where constants q, i = 1,...,7, permits us to carry the SM particles in the first doublet and 
the information about the new physics by the second one. The Yukawa interaction 

U,FC = ifJQ^L^2U,R + i^,Q^L(t>2D,R + h.C. . (6) 
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describes the Flavor Changing (FC) one beyond the SM. Here, the couphngs ^^'^ for the 
charged FC interactions are 

i^h = VCKM , 

^ch = VCKM ) (7) 

and 

where the index " N" in ^'^'^ denotes the word " neutral" . Notice that Hi and H2 are the mass 
eigenstates and respectively, since no mixing occurs between two CP-even neutral bosons 
and in the tree level, for our choice. 

The exclusive t^t^ process is induced by the inclusive h st^t^ decay. Therefore 

we start with the effective Hamiltionian of 6 — st^t^ decay 

- 2C7^ (w^.p.Pij&)r7^r + Cq, (sP,j6)fr + Cq, (sPR6)r75r| , (9) 

with p = pi + P2, the sum of the momenta of r"*" and r~. Note that, Ci,Cg and Ciq, are 
the Wilson coefficients normalized at the scale fi and given in Appendix B. The additional 
Wilson coefficients Cg^ and Cq^ are due to the NHB exchange diagrams (see Appendix B). In 
calculating the Tie//, one first integrates out the heavy degrees of freedom, namely t quark, 
W"^, H"^, Hq, Hi and H2 bosons in the present case and then obtain the effective theory. Here 
H^ denote charged, and H^, Hi ,H2 denote neutral Higgs bosons. Note that Hi and H2 are 
the same as the mass eigenstates and A'^ in the model HI respectively. At this stage the 
QCD corrections are added through matching the full theory with the effective low energy one 
at the high scale /i = mw and evaluating the Wilson coefficients from mw down to the lower 
scale /i ~ O(mfe). 

In the model HI, the neutral Higgs particles bring additional contributions (see eq.(|T9D) 
since the mass of r lepton or related Yukawa coupling ^^^^ enter into the expressions (see 
||14||). Finally the neutral Higgs boson (NHB) contributions are calculated using the on-shell 
renormalization scheme to overcome the logarithmic divergences. Using the renormalization 
condition 

^neutriP ) ~ ^neutriP ) + ^neutr = (10) 
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the counter term T^^^^^ and then the renormahzed vertex function r^g"j^(p^) is obtained. Here 
the phrase neutr denotes the neutral Higgs bosons H^, and and p is the momentum 
transfer. 

For the exclusive decay Bg t^t^ , Tieff is to be taken between vacuum and |-B° > state 
as < 0\Heff\Bg > and this matrix element can be expressed in terms of the decay constant 
Ibs using 



< 0|S75&|5° > 



^fE 



0. 



m 



Bs 



(11) 



Since p = pi +P2, the Cg-^-^ term in eq.(P) gives zero on contraction with the lepton bilinear, C7 



gives zero by eq.([TI|) and the Cio term gets a factor of 2111^- while the remaining Cq^ and Cq^ 
terms get m^^, when taking ^ rnfy + m^. Thus the effective Hamiltonian eq. results in 
the following decay amplitude for Bg 

Gpa 



A 



2^/271 



ts 



Cq, TT 



Tfl 

(Cq, - 2^C,o)fl5T 



(12) 



To calculate the branching ratio we find the square of this amplitude, then sum over the spins 
of the r leptons and finally integrate over the phase space. This straightforward calculation 
gives for the branching ratio of Bg t^t~ 



BR 



647r3 



rnljBsflJ^tbVt 



is I 



\ 



1-4- 



m. 



'1-4 



\Cr 



Q2 



2 Ciol 



.(13) 



3 Discussion 



In the multi-Higgs doublet models, there are many free parameters, such as masses of charged 
and neutral Higgs bosons and the Yukawa couplings. In the present work we study our process 
in the general 2HDM, so called model HI. The Yukawa couplings, which are entries of Yukawa 
matrices can be restricted using the experimental measurements. In our calculations, we neglect 



all Yukawa couplings except ^%tty ^Nbby ^Ntt respecting the CLEO measurement |T5 



BR{B ^ X,7) = (3.15 ± 0.35 ± 0.32) 10' 



(14) 



This section is devoted to and Cn,tt dependencies of BR for the exclusive decay Bg 



T^T , restricting jCy-^'^^l in the region 0.257 < \C^^^\ < 0.439 due to the CLEO measurement. 
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eq.(|T^ (see [!€.] for details). In our numerical calculations, we take the charged Higgs mass 
= 4:00 GeV and the scale fi = mf,. Further, we use the redefinition 



U,D 



and the input values given in Table (|I]). 



V2 



Parameter 


Value 


Tflj- 


i . ( V ) 


Tflc 


1 1^ Vje V J 


rub 


4.8 (GeV) 


rriHO 


150 (GeV) 


rriho 


80 (GeV) 


rriAo 


80 (GeV) 


"em 


129 


A* 


0.04 


mt 


175 (GeV) 


mw 


80.26 (GeV) 


mz 


91.19 (GeV) 


J^QCD 


0.225 (GeV) 


ocs{mz) 


0.117 


sinOw 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 



Fig. m shows 



mfc 



dependence of the BR of the decay under consideration for ^^^^ 



200 GeV^ and the ratio \rtb\ = \jrr^\ < 1- The BR is restricted to the region bounded by solid 

?JV,tt 

lines for C^^^ > or to the small dashed line for Cj^^ < 0. This quantity is sensitive to 



ruh 



and it increases by an amount %60 in the interval 20 < < 80. Besides, the enhancement 



compared to the SM case is predicted as beeing %80. The dependence of the BR for 



mi, 



Ttb > 1 is presented in Fig. |^. For this case, the BR increases considerably even for small 
values of ^^^t-, which is taken 20 GeV, in this calculation. The BR enhances with increasing 

^^T-T, especially for Cj^'^ < case. This figure shows that the BR is strongly sensitive to the 

FD 



parameter ^j^bb ''"tb > 1 and it may get the values four (two) times larger compared to the 

r,N,bb 

.N,rr for l^fel < 1 



ones in the SM for C^^^ < {Cj^^ > 0) even at ^^^^^ = 2 rub. 

represent the dependencies of the BR on the parameter 



Figures 

and rtb > 1 respectively. In rtb < 1 case, the BR increases almost 1.5 times compared to the 
one in the SM for large values of tt = 500 GeV^ (Fig. ^). However, for rtb > 1, this 

enhancement is quite high as shown in Fig 



Even for small values of ^^bb ^ 



N,TT there 



is a possible increase nearly (more than) one order of magnitude compared to the SM case for 
c^^f > {C^" < 0). 

Now we would like to summarize our results: 

• There is a possible enhancement in the BR at the order of magnitude % 150 for \rtb\ < 1 in 
the model III compared to the one in the SM for large values of the model III parameters, 
^N,bb — SOmfe and ^^^^^ = 500 GeV^. The BR is not so much sensitive to the model 
parameters given above. Further, the NHB effects become sizable with increasing values 

• For Ttb > 1, there is a considerable enhancement at the one order of magnitude compared 
to the SM, even for the small values of C^m Cn,tt- ^^i^ case, the BR is larger and 
more sensitive the model parameters for C^^-^ < than the ones for C^^-^ > 0. Note that 
the enhancement for the increasing values of the ^^.^^ is due to the dependence of the 
BR on the NHB effects. 
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Appendix 



A The operator basis 

The operator basis in the 2HDM (model III ) for our process is [l^, [l^, |l^ 



02 = {sLa1^lCLa){cLpl^hL|3), 

03 = {sLali^bLa) J2 i^LfsYgLp), 

q=u,d,s,c,b 

q=u,d,s,c,b 
q=u,d,s,c,b 

Oq = {sLa1t,hLp) iqRfdl^QRa), 
q=u,d,s,c,b 

g 

g 

g 

q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 

Qe = Yf^(^l&^) E i^RQl) 

q=u,d,s,c,b 

Q7 = ^2(''1^"'^r) E 

q=u,d,s,c,b 

Qs = ^2(''l^"^bR) E (4^, A) 

q=u,d,s,c,b 



Q9 = ^2(''L^"'b^R) E ifL^.uq'^R) 

q=u,d,s,c,b 



i2 



Qio = 4^2(-'l^"^bR) Y: (15) 

q=u,d,s,c,b 

where a and (3 are S'?7(3) colour indices and JF^*^ and Q^'^ are the field strength tensors of the 
electromagnetic and strong interactions, respectively. Note that there are also fiipped chirality 
partners of these operators, which can be obtained by interchanging L and R in the basis given 
above in the model III. However, we do not present them here since corresponding Wilson 
coefficients are negligible. 

B The Initial values of the Wilson coefficients. 

For the sake of completeness we also give the initial values of the Wilson coefficients for the 
relevant process. In the SM they are [0 

C^^'imw) = 1, 



■iSM/ \ —x^ + 5xf + 2xt 



4(xt - 1)4 2A{xt - 1] 

,-y>3 I C^,-y>2 



Cf*'(m„.) = —^B(x,) + L:±^C{x,} - D(x,) + |, , 

Cfo^'irnw) = -4^iB{x,)-C{x,)) , 
sm 6w 

C^^imw) = z = l,..,10. (16) 
The initial values for the additional part due to charged Higgs bosons are 

Ctsimw) = 0, 

C^{mw) = Y'F,{yt) + XY F^iyt) , 

C^{mw) = Y^G^) + XYG^ivt) , 

Ci'.imw) = Y'LM, (17) 



where 

Y _ 

rrib \ ' Vtb 



V — ^ f FD j_ FD ^ts 
^ — _ I <,N,bb "T <,N,sb 



Y 



1 

mt 



V* 



N,tc 



(18) 



ts , 



and due to the neutral Higgs bosons are [|I4| 



^A°ffFU \2- 
^Q2\\'iN,tt) . 

^Q2 \^N,tt. 



+ 
+ 
+ 



{{yt - l)(-ei(zA) + {yt - 1)za) + Qi{za) In 2/0, 

0l(^A)l^ yt{xy + ZA)\ 



2ZA- 



yt-i 

_y I UL 



02(za) ' 63(2^)- 



+ 2(1 + In 



+ 



za 



QlizA) 



-e3(zA)+(x-y)(xt-yt)zA 



m 



w 



Aza ^ 2{x{xt + yt) - 2yt)zA 2{xt{x - 1) + (x + l)yt)zA 



m 



AO 



Q2{za) 



-Qz{za) 



Qz{.za) + {,x- y){xt - yt)zA 



{yt - Xt)zA + Xtyt{2zA - 1) (4a;? - 7yt - 4x?(2 + yt) + xt{5 + 8yt + f-)) Inxt 



{xt - l){yt - l)zA 
{xt{^-l)-yt)lnyt 



{xt - lY{xt - yt) 



{xt - yt){yt - ^)'^ 



+ 41n 



\Q2{za) 



za 



^N,bb) 



9 ^N,TT^N,bb f '^QsizA) - xt{{x - 2)yyt - Xt{yt - za)) 



.4" 



QsizA) 



2 In 



QsizA) 



Za 



{yt - 2xt{yt + 1) + x?(^ + 1) hix^ ^ (x,(l - 2y t) + 2yt{yt - 1) + x?(f^ - 1)) ln7/t\ 



{xt - l){xt - yt) 



iVt - l)ixt-yt) 



9 {^N,u) rnbuiryt / e^za) + (yt-i) 



-l+4yt+yf{2lnyt-3) 



H'^llcU \2\ 



N,tt) 



2bQTi'^m'\jom^rXt \ 



cos^ 6w 



2zh{-Qa{zh){1 - 2x)xt + 2ei{zH)x) {yt - 1)(1 + + yt{xt - 3)) + 2yt{yt - Xt) Inyt 



Qi{zh)Q^{zh) 
9^^N,tt^§,bb^r (yt{2-xt)-xt 



NM) 



GAir'^m'jjomt 
zh 



yt-i 



+ 2xt In 



{yt - 1)3 

rei(zH)i {xt{l- hyt) + 2yj{l + Xt)) \nyt 



Zh 



{yt - If 



Q-^{^zh)Qa{zh)^ - x^xt<c)-,{yt - 1) + x{xt{-<c)i{zH) + 65(1 + X - ?/))(?/t - 1) 

+ 2yt{-Q^ + 2y2/i)) + {-xtQ^{yt - 1) - 2(1 + 2/(2/^ - l))yt)zH) 
yt{{yt - l){-Qi{zH) + zh{i - yt)) + Q>i{zH)yt Iny*)- 



+ 



cos^ 6'vy04(2;/f)(yt - 1)^ 

/ 4,(x-l) I xt(xt(4-xt)-3+2a;f (3;t-2) Inat) 



67 



+ 



bl2Ti'^m\jQmy^ 



cos^ 6 



w 



4(x(4 + 2xt - ff ) - 2) 4(x4( 



3 + 67) - 4x) + 207Xi In 67 



eg 2:^97 

2(2 - 12xj + 21x^2 - 12x^3 + a;4 + (2 - 4xt - 2x2 + 6x?) Inx*) 



-4xt(l + 21n 



6AQi{zh)ml„'mt7r'^{-l + ytY 
2{2x -l){yt- If Zh)-2Qi{zh)yt\^yt) 



((-i + yt)(e>i(-2fe)(yt + i) 



nholItU ^2^ _ ^ cD cD (?U ^2 ( ^ijZh) + Vtjxy - Zh) , (1 - j/t + (- 1 + 2|/t) 111 |/t 



2 In 



). 



^/.o -u _ g^CN,rAN,uXt I xt(8 - 9yt) - (y* - 2) + yt{hyt - 4) + (-4 + 2yt + y|) 



1287r^m^j(,mj 



OsCxi - 1) 



ytXt{2 - xt - yt) 4zh{-l + x{2 + Xt)) - 2xyxt , 22;^(-2|/t + a;(a;t + y*)) 



^ 22:^(0:^(0: (a: +1)) ^ 



^ ©3(^0 
4(-l + x)xxtytZh 



+ 



-©3(^/1) + (a: - i/)(o:t - yt)Zh {Qi{zh)xt - x{xt - yt)zh){Qi{zh)xt - y{xt - yt)zh) 
. {yt-lf{xl{?>-l^yt) + 7y^,-7xm{'^ + yt) + ^xl{l + Ayt + 2yl)) 

ytXt{-l + ytfj-Oe + 4(a;t - yt)yt) \ Inxt 
ZH ^ ©i 

{x, - lf{-lQx,yt{yt - 1) + xl{2y, - 1) + ^^(4^, - 5) - ^©6)lny, _ ^ ^e2(^.) 



©^ 



Zh 



). 



g'ilrr^^ZH ( (xf (t + 1) + - 2o:,(y, + 1)) In X, 



where 



QAn'^m^rXt 



{xt - l){xt-yt) 



+ 



1 2;? [^^ - 1) + 2t/,(t/f - 1) + xt(l - 2yt)) Inyt 



{yt - l){xt - yt) 



-2 In 



©3(-2^) 



Zh 



yt{xt + 2y{x-l))-Zh{xt-2yt{y-l) + 2yf) + x{yyt + 2zh{yt-l))^ 



©2(0;) 
©3(0;) 

©4(0;) 

©5 
©6 
©7 



©8 = 



((1 - 2/ + - x{yyt + a;(l - yt)) 

©i(^,yt ^ 2:4) 

(^^((l -y)+ y)yti^ - xxt{yyt + u;{-l + yt)) 

(2/(1 - yt) + yt)^ - xiyyt + uj{-l + yt)) 

(-1 +xt){xt -yt){-l + yt) 

{-I + y){y{-l + yt)yt) 

{xt + Z/(l - + o:(2:fe - Xt{y + 2:^)) 

(1 - y{l - xt))zh - x{xt{y - Zh) + Zh) 

XtZh 



(20) 
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and 



yt 



m; 



mH± 



mr 



Zh 



m 



Zh 



HO 



m 



za 



m 



(21) 



AO 



The explicit forms of the functions Fu2){yt), Gi(2){yt), Hiiyt) and Li{yt) in eg. ([T7|) are given as 



Fiivt) 

G2{yt) 
Hiiyt) 



yti7 -5yt 



iyt. 



+ 



72(y,-l)3 I2{yt-IY 
ytihyt - 3) yt{-^yt + 2) 



12(1/* 



+ 



6(?/t 



l/i(-l/,^ + byt + 2) 
24(|/i - 1)3 

yt{yt - 3) , 



-yt 



^{yt 



+ 



4(|/t - 1)2 2{yt - 1 



Asiri^Ow xyt 



yt 



t 



1 



1 



yt - 1 (yt - 1)^ 



108(i/t - 1)= 



Li{yt 



xyt 



yt 



+ 



iKyt - 1] 
1 



(22) 



Finally, the initial values of the coefficients in the model III are 



a 



2HDAI^ 

2HDM I 
h 



Hi 



dx 



l-x 



dyiGg^ii^N., 



hO i?U 



ihO i?D 



ri2HDM 

^Q2 



{mw) 



l-x 



dx I dy (C4°((e^,j3) + cg^^^l,,f) + Cg{eN,t) + G^:{eN,b)) 



ri2HDM(^ \ 



G, 



2HDM I 



mw] 



G, 



2HDM , 



mw] 





rrib 


iGl^.^'^'irnw) 


1 r^2HDM 
+ *-^Q2 


mr 


sin^ 9w 




rrib 


{Glg'"\rnw) 


ri2HDM 


mr 


sin^ 9w 


= 0, 


i = 5, . 


.,10. 





(23) 



Here, we present Cq^ and in terms of the Feynmann parameters x and y since the inte- 
grated results are extremely large. Using these initial values, we can calculate the coefficients 
(j2HDM^^^ and CQ^'^*^(/i) at any lower scale in the effective theory with five quarks, namely 



u, c, (i, s, b similar to the SM case |]13|, [T^, |T9|, |22|. For completeness, in the following we give 
the explicit expressions for G^"{fi) and Gl^^{n). 



G'r"{f^) 



Cr^''(/i) + {Gr^'^'if^) + iVc Ce^™'(/x)) 
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+ g„ i'-^ c^r'Hf^) + c^r^'if^)) , (24) 

nib rrib 



where the LO QCD corrected Wilson coefficient is given by 

^LO,2HDM^^~^ = 7716/23^2™/^^^) ^ (8/3)(r7l^/23 _ ^16/23)^2/fDM(^^) 
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+ C2™^(mH.)EM"' . (25) 



i=l 



and ?7 = as{myY) / Cis{lj), hi and Oj are the numbers which appear during the evaluation 
The Wilson coefficient C^^^ {^) is : 

+ h{z, s) (3Ci(/x) + C2(/x) + SCsifi) + C.ifi) + 3C,{fi) + CM) 

- s) (4C3(^) + 4C4(/i) + 3C,{fi) + C6(/i)) (26) 

- ^/i(0, s) (Csifi) + SCifi)) + ^ (SCM + Cifi) + 3C,ifi) + CM) ■ 

Here the functions h{u, s) are given by 

hf \ 8 mfe 8 8 4 , . 

n{u,s) = — m lnu-\ — x (27) 

^ ' ^ 9 ;u 9 27 9 ^ ' 



/l~x-l 



- in] , for a; = — < 1 



9 I 2arctan^i=^, for x = — > 1, 



with M = 

mi, 



Finally, the Wilson coefficient Cio(/i) is the same as Cio{mw) and Cq-^^fi), Cq^i^fi) are given 
by El 

C^Q.(/i) = r''/''C'Q,(mH/) , ^ = 1,2. (29) 
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Figure 1: BR as a function of ^N,bb/''^b for ^^.^^ = 200 GeV in case of the ratio \rtb\ < 1- 
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Figure 2: Si? as a function of ^N,bb/''^b for C^,tt — 2,0 GeV in case of the ratio rtb > 1- 
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